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Abstract 



Data taken by the DELPHI experiment at centre-of-mass energies of 183 GeV 
and 189 GeV with a total integrated luminosity of 212 pb~^ have been used 
to search for the supersymmetric partners of the electrons, muons, and taus 
in the context of the Minimal Supersymmetric Standard Model (MSSM). 
The decay topologies searched for were the direct decay — *• -^x?)) Pro- 
ducing acoplanar lepton pairs plus missing energy, and the cascade decay 
{^ — > ^X2 ~^ ^lXi)i producing acoplanar lepton and photon pairs plus missing 
energy. The observed number of events is in agreement with Standard Model 
predictions. The 95% CL excluded mass limits for selectrons, smuons and staus 
are mg < 87 GeV/c^ , < 80 GeV/c^ and nif < 75 GeV/c^, respectively, for 
values of /x=-200 GeV/c^ and tan/3=1.5. 
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1 Introduction 

During the 1997 and 1998 data taking period, the LEP accelerator operated at centre- 
of-mass energies of 183 GeV and 189 GeV respectively. This allowed an extension of the 
searches for scalar partners of electrons, muons, and taus, predicted by supersymmetric 
models, over the limits on the production of these particles obtained from data previously 
taken at centre-of-mass energies of 130-172 GeV This paper reports on a search for 
these particles using the 212 pb~^ of data taken by DELPHI during 1997 and 1998. Similar 
searches have been performed by other collaborations 0. 

For a realistic experimental search one has to make some well motivated assumptions. 
In this analysis, the model assumed is the Minimal Supersymmetric Standard Model 
(MSSM) 0]. In the case that the MSSM is locally invariant (often referred to as minimal 
supergravity) , the number of free parameters set at the unification scale (the scale at 
which gauge couplings unify) can be reduced to fiveQ: 

mo, mi/2. A, B, /i. 

These are respectively, the universal scalar and gaugino masses, the universal trilinear 
and bilinear^ scalar couplings, and the Higgs doublets mass mixing parameter. 

In this analysis R-parityQ conservation is also assumed, which leads to three important 
phenomenological consequences 0. Firstly, the lightest supersymmetric particle (LSP) 
must be absolutely stable. If the LSP is electrically neutral, as favoured by cosmolog- 
ical constraints, it interacts only weakly with ordinary matter, so escaping detection. 
Secondly, each supersymmetric particle (sparticle) other than the LSP must eventually 
decay into a state which contains an odd number of LSPs, typically just one. Finally, 
R-parity conservation implies that collider experiments could only produce sparticles in 
even numbers. 

Consequently, sleptons could be pair produced at LEP via e+e~ annihilation into Z°/7 
(Figure l.a). In addition, selectrons can be produced from t-channel neutralino exchange, 
which introduces a direct dependence on the SUSY parameters and the possibility of left 
and right handed final states even without mixing via the mass matrix (Figure l.b). 

In a large fraction of the SUSY parameter space the dominant decay of the sleptons {£) 
is to the corresponding lepton flavour plus the lightest neutralino (x?) (Figure 2. a), pre- 
sumed from the MSSM mass spectrum to be the LSP. The neutralino will escape unde- 
tected, hence the topology will be characterised by acoplanar lepton pairs together with 
missing energy. In most of the analyses described in this paper we search specifically for 
such a signature. 

The search can be extended by looking for topologies other than acoplanar lepton pairs. 
For certain values of SUSY parameters it is possible for the second lightest neutralino X2 
to be lighter than the sleptons. If this is the case the slepton can also decay via a cascade 
to a Xi5 with a possible decay chain i ix2 ^ ^IXi (Figure 2.b). In order to investigate 
this channel we have searched DELPHI data for events containing acoplanar lepton and 
photon pairs plus missing energy. 



^See H and references therein for further information on the actual supersymmetry breaking mechanism, and 

motivation for tne assumptions made. 

^ Using renormalisation group evolution, the bilinear term is expressed in the low energy MSSM as tan/3, the ratio of 
the vacuum expectation values of the two Higgs doublets. 

^R-parity is a quantum number, defined as R = jgj^ with B, L, and S respectively the baryon number, 

the lepton number and the spin of the particle. Non-supersymmetric particles, including the Higgs scalars axe H-even, 
whilst the supersymmetric particles are all K-odd. 
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2 Detector description 

The DELPHI detector and its performance have been described in detail elsewhere [|^ ^j; 
in the following we present only a brief description of the components relevant to the 
analyses presented here. 

A system of cylindrical tracking chambers coupled with a 1.2 T uniform solenoidal 
magnetic field, directed along the beam axis, enables the reconstruction of charged par- 
ticle tracks. The Vertex Detector (VD) consists of three cylindrical layers of silicon 
detectors, at radii 6.3 cm, 9.0 cm and 11.0 cm. The vertex tracking is aided in the 
forward regions by mini-strips and pixel detectors making up the Very Forward Tracker 
(VFT) with an angular acceptance between 10° and 25°. The Inner Detector (ID) 
is a cylindrical drift chamber (inner radius 12 cm and outer radius 22 cm). The Time 
Projection Chamber (TPC), the principal tracking device of DELPHI, is a cylinder of 30 
cm inner radius, 122 cm outer radius and a length of 2.7 m. Each end-plate is divided 
into 6 sectors, with 192 sense wires used for the dE/dx measurement and 16 circular 
pad rows used for 3 dimensional space-point reconstruction. The Outer Detector (OD) is 
composed of 24 planks each with 5 layers x 32 columns of drift tubes. The tubes, situated 
at radii between 196 cm and 207 cm from the beam axis, improve the precision of the 
momenta of the charged particles measured by the TPC. In addition to the barrel track- 
ing, two planes of drift chambers. Forward Chambers A (FCA) and B (FCB), aligned 
perpendicular to the beam axis, allowed tracking in the endcap of the detector, giving a 
polar coverage down to 11° and 169° with respect to the e~ beam direction. 

The electromagnetic calorimetery consists of the High density Projection Chamber 
(HPC), covering the barrel region of polar angle 9 in the range 43° < 9 < 137°, the 
Forward ElectroMagnetic Calorimeter (FEMC), consisting of 9064 Cherenkov lead glass 
blocks covering 11° < 9 < 36° and 144° < 9 < 169°, and the STIC (Scintillator Tile 
Calorimeter), extending the coverage down to 1.66° from the beam axis in either direction. 
The 40° taggers are a series of single layer scintillator lead counters used to veto photons 
and electrons that would otherwise have been missed in the region between the HPC and 
FEMC. 

The hadron calorimeter (HCAL) covers 98% of the solid angle. Muons with momenta 
above 2 GeV traverse the HCAL and are recorded in a set of muon drift chambers; 
the MUon Barrel (MUB) chambers, MUon Forward (MUF) chambers and the Surround 
Muon Chambers (SMC). 

The identification of muons is provided primarily by the algorithm described in 0, 
which relies on the association of charged particles to signals in the barrel and forward 
muon chambers. In order to reduce contamination from cosmic ray particles, the impact 
parameter with respect to the beam crossing point was required to be less than 1.5 mm 
in the R — cf) plane. 

Electrons are identified as charged particle tracks with an energy deposit above 3 GeV 
in the electromagnetic calorimeter and with the ratio of the electromagnetic calorimeter 
energy to the track momentum from the track above 0.3c. In addition, the shape of the 
shower profile in the HPC and the dE/dx measurement in the TPC were also considered. 
Forward electrons are distinguished from gamma conversions by requiring hits in the 
VFT. 

A charged particle is identified as a pion if the energy deposited in the HCAL is greater 
than 2 GeV, greater than the energy deposited in the electromagnetic calorimeters and 
it does not produce hits in the muon chambers. 
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3 Data samples and event generators 



The total integrated luminosity accumulated by the DELPHI experiment over the two 
years analysed was 212 pb^^. This included 54 pb^^ of data collected at a centre-of-mass 
energy of 183 GeV and 158 pb"^ collected at 189 GeV. 

Several programs were used to simulate Standard Model (SM) and SUSY (signal) 
events in order to estimate background contamination and signal efficiencies. 

All the models used JETSET 7 . 4 |jl^ for quark fragmentation with parameters tuned 



to represent DELPHI data []TT|. The program SUSYGEN [|T^] was used to generate slepton 

events and to calculate cross-sections and branching ratios. The generator EXCALIBUR |jl3| 
was used to model all four-fermion events, which includes the coherent interference of all 
diagrams leading to a given final state. For a cross check, PYTHIA was used to 
generate samples of WW, ZZ, Wez/ and Zee events. The processes e+e~ — »Z°/7 — > qq(7) 
were simulated by PYTHIA, whilst the two-fermion backgrounds e"'"e~ —^T?j^ fi^fi^{'-f) 



and r+r (7) events were produced by KORALZ |T1]. The generators BABAMC [|T3| and 



BHWIDE [O] were used to simulate Bhabha scattering. Two-photon interactions leading 
to hadronic final states were simulated using TWOGAM [|1^ and BDKRC ||18| for the Quark 
Parton Model contribution. BDK [|T^ was used for final states with electrons only, whilst 
final states with muons or taus were simulated using BDKRC. 

Generated signal and background events were passed through a detailed detector re- 
sponse simulation (DELSIM) @] and processed with the same reconstruction and analysis 
programs as the data. The number of background events simulated was several times 
larger than the number expected in the data. 



4 Search for select rons and smuons (i ^Xi) 

The analysis was performed in two stages. Firstly a loose pre-selection was used to 
obtain a sample containing events with two oppositely charged tracks. At this stage, 
various distributions of the real data were compared with distributions from simulated 
SM events. 

After this stage a tighter selection was applied. Tuned to both simulated background 
events and signal events, selections were made in order to reduce the expected SM back- 
ground whilst keeping a reasonable efficiency for the signal over a wide range of the i-Xi 
mass combinations. 



4.1 Search for selectrons 

To search for selectrons, the general topology required was two acoplanar electrons 
and missing energy. The preliminary event selection kept all candidates with exactly two 
well reconstructed oppositely charged particles with momentum above 1 GeV/c. One of 
the two charged particles was required to be identified as an electron, rejecting events if 
the other was identified as a muon. At this stage in the analysis, the selection consisted 
mainly of Bhabha and two-photon events, with satisfactory agreement observed between 
data and simulated background (Figure 

To further reduce the SM backgrounds, tighter cuts were applied. As two-photon 
events are predominantly at low polar angles and with low momentum it was required 
that the visible energy be greater than 15 GeV and that the energy deposited in the low 
angle STIC calorimeter be less than 4 GeV. As a further constraint, the invariant mass 
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of the two tracks was required to be greater than 4.5 GeV/c^, and the total transverse 
momentum with respect to the beam axis was required to be greater than 5 GeV/c. 

To reduce the number of Bhabha events an upper hmit on the visible energy of 100 GeV 
was imposed, whilst also requiring that the neutral energy not associated to the charged 
tracks be less than 30 GeV. Events were also rejected if there were more than four neutral 
clusters in total, each with energy above 0.5 GeV. Bhabha events are coplanar with a 
large opening angle, hence it was necessary that the acoplanarity and acolinearity be 
greater than 15°. 

Four-fermion events were reduced by the constraints described above, in particular the 
constraint on the visible energy. 

Contraints were also imposed on the momenta of the two tracks, requiring that both 
tracks had momenta above 2 GeV/c. It was further required that the missing momentum 
vector pointed to an active region of the detector. 

The efficiency for the signal detection depends on the masses of e and Xi- The typical 
signal efficiency is ~ 50 %. 

After this selection a total of 56 candidates were found in the 212 pb~^ of data analysed, 
compared to 51.2 ± 1.5 predicted from SM processes. Details are given in Table |l|. 

4.2 Search for smuons 

As a pre-selection, exactly two well reconstructed oppositely charged particles with 
momenta above 1 GeV/c were required. At least one of the particles had to be identified 
as a muon. It was further required that neither track be identified as an electron. The 
pre-selection sample consisted mainly of two-photon events, and good agreement between 
real data and simulated background was observed (Figure |^). 

To further reduce SM backgrounds tighter cuts were applied. It was seen that using 
a sequential cut analysis, the dominant background after a tighter selection was W-pair 
events. These events become increasingly important in regions of high slepton mass and 
high Am (m^-m~o), where the signal events become virtually indistinguishable from the 
W-pair background. In these regions of SUSY mass space the cross-section for smuon 
production is lowQ, and hence using sequential cuts to remove this background has a 
severe effect on the signal sensitivity. 

Consequently for the 189 GeV data analysis, a different approach was adopted to 
the 183 GeV analysis, applying a selection procedure which depended on the (/i, Xi) 
mass difference. For a mass difference. Am, less than 35 GeV/c^ where two-photon 
backgrounds are important, an analysis based on sequential cuts was performed. For 
the data taken at 183 GeV this approach was used for the full SUSY mass spectrum. 
However, in the 189 GeV analysis, for regions of Am greater than 35 GeV/c^, where 
the W-pair backgrounds are kinematically favoured, a probabilistic analysis based on the 
likelihood of an event being compatible with W-pair production was used. 

For regions of Am less than 35 GeV/c^, to remove the two-photon events, the visible 
energy was required to be greater than 10 GeV. Also the energy in the STIC had to be 
less than 1 GeV. As a further constraint it was necessary for the invariant mass of the 
lepton pair to be greater than 4.5 GeV/c^. 

To remove e+e~ ll events in this Am region an upper limit on the 

visible energy of 120 GeV was imposed whilst also requiring the unassociated neutral 
energy to be less than 10 GeV, with no more than two neutral clusters. This background 

'*In the selectron scenario the t-channel contribution can enhance selectron production for low neutralino masses, hence 
increasing signal sensitivity in the mass regions dominated by W-pair backgrounds. 
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was further suppressed by accepting only events in which the opening angle between the 
tracks was less than 165° and the acoplanarity was greater than 15°. 

To reduce W-pair contamination in this low Am region, at a small cost in signal 
efficiency, events were rejected if the positively charged muon was within 40° of the e"*" 
beam direction, or the negatively charged muon was within 40° of the e~ beam direction. 

For the selection of events kinematically allowed in regions of Am greater than 
35 GeV/c^, a discriminating variable was constructed for the events in the 189 GeV 
data using the probability density functions (p.d.f 's) of W-pair event variables after the 
pre-selection stage. The following variables were chosen due to their high discriminating 
power between signal and four-fermion events and their relatively low correlations: 

• Product of lepton polar angle and charge (Q cos6'); 

• Neutral energy; 

• Opening angle between the leptons; 

• Acoplanarity; 

• Missing energy; 

• Missing transverse momentum. 

In addition, these variables have excellent agreement between real data and simulated 
background. 

The discriminating function is shown in Figure 5 for data and simulated background 
(which is predominantly 2-photon at this stage). Also shown is the comparison of the 
discriminant variable for a sample of 4-fermion events and a sample of SUSY signal with 
a Am value close to the W mass (80 GeV/c^). 

Two-photon and di-muon events were removed using the same cuts as in the low mass 
window. W-pair contamination was reduced by cutting on the discriminating function 
such that signal to background was maximised. It was further required that the missing 
momentum vector pointed towards active components of the DELPHI detector. 

The efficiency for the signal detection depends on the masses of ft and Xv The cuts 
used to remove the SM background resulted in typical efficiencies of ~ 50 % for the 
regions of low Am, and ~ 35 % for the regions of high Am. 

Table ^ summarises the number of accepted events in the data together with the 
predicted number of events from background sources. In the data collected at 189 GeV, 
for the regions of Am< 35 GeV/c^, 17 candidates passed the tight selection, consistent 
with a background prediction of 17.5 ± 0.3 events. For the regions of Am > 35 GeV/c^, 
7 candidates remained compared with a background prediction of 9.2 ± 0.2 events. In 
the analysis of data taken at a centre-of-mass energy of 183 GeV, 5 candidates remained 
with an expectation of 6.1 ± 0.6 events from SM processes. 



5 Search for staus 

The off-diagonal terms of the slepton mass matrix are proportional to the mass of the 
corresponding SM partners. Important effects caused by these terms must be considered 
when searching for staus. For a certain mixing between right and left handed staus, 
the low-mass stau eigenstate (fi) can become an important candidate for the lightest 
charged supersymmetric particle. Another consequence of the possible mixing of staus 
is the change of the coupling to the Z°, and consequently the production cross- sect ion, 
with the mixing angle. 
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5.1 Search for heavy staus 

The characteristic signature of the production of pairs of heavy staus is the detection 
of a r+r~ pair with large acoplanarity and missing energy. Due to the scalar nature of 
the stau, the visible system will tend to be at large angles to the beam. 

Among the SM background processes to this signal are s-channel production of tau 
pairs, in particular if they arise from a radiative return to the Z°, with the ISR photon 
escaping detection, and four-fermion events where the final state contains two taus as the 
only visible particles. Finally, two photon interactions with 77 — >■ t~^t~ contribute in the 
case of staus close in mass to the LSP. 

To select events with two taus, well reconstructed charged and neutral particles were 
first collected into clusters of total invariant mass below 5.5 GeV/c^. Events with exactly 
two particle clusters (possibly accompanied by isolated neutral particles) were considered 
further if there were no more than 6 charged tracks in the event and these gave a total 
charge of or ± 1. At least two tracks were required to have momentum above 1 GeV/c 
with one greater than 4 GeV/c. The distributions of data and simulated SM events agree 
well at this stage in the analysis, as can be seen in Figure ^ 

To ensure that the selected events had the high acoplanarity typical for the signal, the 
acoplanarity angle was required to be above 10° (11° for the 183 GeV sample). 

Selecting events at high angles to the beam was done by demanding that at least two 
charged particles with momentum above 1 GeV/c were observed above 30° to the beam 
axis. Also, the direction of the vectorial sum of momenta should be contained in the 
barrel region at an angle greater than 37° to the beam (30° in the 183 GeV sample). 
To reduce the background from radiative returns to the Z°, none of the clusters were 
allowed to have a total momentum {p'^^'^) above 67 GeV/c (60 GeV/c in the 183 GeV 
sample), the energy of isolated photons had to be below 20 GeV and there should be no 
signal in any 40° tagger. Furthermore, the value of the reduced centre-of-mass energy 
{Vs') was estimated using the angles of the jets, and assuming that a photon was lost 
in the beam-pipe (the triangle rule). This value should not fall in the interval 90 to 94 
GeV (no such cut was made for the 183 GeV sample). In addition, there had to be no 
calorimetric energy below 30° in polar angle. This last cut was also very effective against 
the background from two-photon events, as it removed all such events where either of the 
initial e"'"e~ were deflected into the forward calorimeters. The two-photon background was 
also reduced by discriminating against events where the two clusters were close together: 
the acoplanarity angle should not exceed 170° (176° for the 183 GeV sample). 

In order to further suppress the background from e"'"e~ — >■ Z°/7 — > r+r" events with 
r-decays highly asymmetric in visible momentum, the square of the transverse momentum 
with respect to the thrust axis (6) had to be above 0.9 (GeV/c)^. This was the case for the 
189 GeV sample; the condition was more complex in the 183 GeV sample (see below). 

At this stage of the analysis, the background was dominated by the W-pair background, 
and in order to further suppress it, the events were analysed under the assumption that 
they were indeed W-pair events. The 9 angle of the positive W {Oiy+) was reconstructed 
by assuming that the direction of the taus was indentical to that of the jets and applying 
an unsmearing procedure (derived from simulated W-pair events) to estimate the true 
momentum of the taus. The final estimate of 6w+ was then given by the average of the 
polar angle of the positive jet, the complement of the polar angle of the negative jet, and 
the two approximate solutions to the equation determining the W angle. As the signal 
is isotropic, whilst the W production is enhanced in the forward-backward direction and 
concentrated at high values of the higher of the two jet momenta (p^^), it was required 
that the observed values of 9w+ (in radians) and p^^x (i^ GeV/c ) were below the higher 
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of the two lines 6w+ = 1-5 and 6^+ = -0.05 Pmlx + 3-7 in the p^n^ — Ow+ plane. This 
cut was applied to the 189 GeV sample; in the 183 GeV sample, the simpler cut Ow+ < 
2.5 radians was used. 

This selection was supplemented by cuts that depended on the region of the (m^ ,771^0) 
plane considered, which were tuned to remove the corresponding backgrounds for that 
region. The sub-division was different in the two samples, as were the cuts depending 
on the region. For the 189 GeV sample, two regions were defined: Am = irif-m^o 
below or above 20 GeV/c^. In the region of low mass difference, where the remaining 
two-photon background was concentrated, it was required that the missing transverse 
momentum (p™'^'^) was greater than 5.4 GeV/c whilst, in the region of high Am, the 
p™*^*^ was required to be above 8 GeV/c. In the lower (higher) Am region, it was also 
required that the highest momentum of any identified lepton in the event was less than 
20 GeV/c (22 GeV/c), in order to further suppress the remaining W-pair background. 

For the 183 GeV sample three regions of Am were considered: less than 22 GeV/c^ 
22 — 50 GeV/c^ and more than 50 GeV/c^. In these three regions, 5 was required to 
exceed 0.4, 1.0 and 0.4 (GeV/c)^ respectively, and p™^^ to exceed 5.5, 6, and 6 GeV/c, 
respectively. The momentum of any identified lepton in the event should be less than 
30 GeV/c (independent of Am). 

Table ^ summarises the number of accepted events in the data for the different se- 
lections together with the expected numbers of events from the different background 
channels. In the 212 pb~^ data sample analysed, 16 candidates were found, with a back- 
ground estimation of 18.1 ± 0.8 from SM processes. The signal detection efficiency was 
of the order of 20% for the 189 GeV sample and 30% for the 183 GeV sample. 

5.2 Search for light staus without coupUng to the Z° 

To a large extent a light stau can be excluded using the agreement of the decay width 
of the resonance with the SM prediction, as observed at LEPl |2^. The corresponding 
cross-section limit of 150 pb for non-standard processes at a/s = Mz excludes a tr below 
25 GeV/c^. However, at the stau mixing angle giving the minimum cross-section, the 
coupling to the TP vanishes and no exclusion is possible using this method. The high 
mass analysis described in the previous section loses its efficiency for stau masses below 
20 GeV/c^. This is mainly due to the fact that the stau-pairs are highly boosted at 
such a low mass, so that they fail the acoplanarity cut. Therefore a specific search 
was required for nif in the range from m,- to 27 GeV/c^ and this mixing angle. After 
selecting two tau events as described in the previous section, it was required that there 
was no identified isolated photon in the event. To reduce the background from Bhabha 
scattering, we demanded that the acolinearity angle was above 0.4° and p'^^^x below 
70 GeV/c. Furthermore, the missing transverse momentum in the event was required 
to be above 6 GeV/c, the angle of the most energetic track in each hemisphere of the 
detector to be above 50° to the beam, and the direction of the vectorial sum of momenta 
should be above a polar angle of 40°. To further reduce the background from radiative 
return to the Z°, the value of \fs' should not be between 82 to 102 GeV 0. 

For TTLf above 15 GeV/c^, the cut on acoplanarity was more restrictive: it was de- 
manded to be above 4°. 

With these cuts, a total of 122 events were selected in the mass region below 15 GeV/c^ 
and 50 events in the higher region. The SM background was estimated to 150.1 ± 2.3 and 
55.7 ± 1.7 in the two regions (Table The efficiencies were practically flat for both mass 

^It should be noted that as a consequence of the vanishing couphng to the Z", radiative return docs not occur in t 
production at the mixing angle yielding the minimal cross-section, even at nif below Mz/'i. 
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ranges. In the first case, the efficiency was around 25%, decreasing to 10% for Am = 171^-. 
In the second case it was around 15%, decreasing to < 2 % for Am< 2 GeV/c^. 

6 Search for cascade decays 

In order to extend the slepton search, topologies from cascade decays of the i have 
been considered. There are regions of the SUSY parameter space where the sleptons may 
also decay into the X2 pl^s the corresponding lepton (i (^xV)- 

For smuons it presents some advantages as the dependence of the smuon mass limit 
on the SUSY parameter /i is considerably reduced. The X2 decay to Xi7; with the 
x\ being the LSP which escapes undetected. The topology for these events is acoplanar 
lepton pairs and two photons plus missing energy. The main advantage of searching for 
this type of event is that the experimental signature is very clean with very small SM 
backgrounds since the emission of photons requires higher orders and extra uem factors. 
However the cascade decay may be suppressed as the branching ratio BR{x2 ~^ Xil) 
may be small. The BR{x2 X1I) is close to 1 when both neutralinos have a similar 
mass, but in this near degenerate case the outgoing photon has low energy, so making 
detection difficult. 

The search for these events was done in a two step procedure. First, samples of 6677 
and /x/i77 events were selected following loose cuts. Before applying a tighter selection a 
likelihood function was defined for the main background channels (one for tagging Bhabha 
events in the 6677 sample and the second one for tagging the e~^e~ — >■ /i"*"/!" entering the 
/^/i77 samples) contributing to the sample. Then the tight selection is combined with the 
likelihood in order to remove those events compatible with the SM processes. 

6.1 ee77 selection 

A sample of 6677 events was selected requiring only two charged tracks reconstructed 
with momentum greater than 3 GeV/c and at least two photons reconstructed with 
energy above 1 GeV. If several photons were selected, only the most energetic two were 
retained. Events were rejected if either of the charged tracks was consistent with positive 
pion or muon identification (see section 2). An acceptance cut demanding that all 4 
particles should lie in the 10° < 6 < 170° region was applied in order to remove most of 
the two-photon interactions. The same cut was applied to the missing momentum vector, 
since this can be close to the beam axis for events coming from a radiative return to the 
Ti'^. Finally the acolinearity of the electron and photon pairs had to be larger than 3°. 
After this stage in the analysis, reasonable agreement was observed between real data 
and simulated SM processes contributing to the sample (Figure |^). 

The Bhabha likelihood was built according to the probability density functions (p.d.f 's) 
of the visible energy, invariant mass between the electron-photon pair with smallest open- 
ing angle, the invariant mass of the electron pair and most energetic photon, and the angle 
between the missing momentum and the closest electron and photon. 

A sample with most SM events removed was selected by demanding that the events 
satisfying the loose 6677 selection comply with the following cuts: acolinearity and acopla- 
narity of the electron and photon pairs above 6° and a Bhabha probability less than 3%. 
According to this selection, 4.2 ± 0.9 events were expected from the SM and 5 seen 
(Table D. 
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The efficiency for signal detection depended on the lepton and photon energies. An 
efficiency map was computed for a range of points of the SUSY parameter space. Typical 
efficiencies for detection were of the order of 40%. 

6.2 AtAt77 selection 

The selection of a sample of /i/i77 events proceeded in a similar manner to that of the 
ee77 events. It was required that at least one of the charged tracks had to be identified as 
a muon, and neither of them as an electron. The rest of the kinematic cuts were the same 
as mentioned in the 6677 selection. Events with more than two photons were further 
considered if the energy of the extra photon(s) was below 10 GeV. 

The e^e^ f^^fJ'" likelihood was built from the p.d.f.s of the visible energy, momen- 
tum of the leading muon, invariant mass of the muon pair with the most energetic photon, 
the angle between the missing momentum and the closest muon, and the opening angle 
of the muon pair. 

The tight fifi'j'j selection consisted basically of the same topological cuts as the tight 
ee77 selection. In addition to these cuts, it was required that the likelihood for an event 
being consistent with e'^e'^ ji^ ji~ be less than 5%. With these cuts 2.9 ± 0.7 events 
were expected from the SM and 3 seen in the data (Table P). The cuts used to remove 
SM background resulted in typical efficiencies for signal detection of 45%. 

7 Results 

Limits on slepton masses can be derived using several different assumptions. Scalar 
mass unification suggests lower masses and cross-sections for the partners of right handed 
fermions. Hence we have assumed that only right handed selectrons {qr) and smuons {jiR) 
are produced, leading to conservative mass limits. 

For third generation sfermions, Yukawa couplings can be large, leading to an appre- 
ciable mixing between the pure weak hypercharge states. The production cross-section 
depends on this mixing, due to the variation in strength of the coupling to the TP com- 
ponent of the weak current, and has a minimum at a mixing angle of 42°. Consequently, 
the results for the stau analyses are presented under these two assumptions; right handed 
stau production (t/j), and minimal mixing stau production (fmm)- 



8 Exclusion limits 

Exclusion limits for slepton pair production were obtained, taking into account 
the signal efficiencies for each i — Xi mass point, the cross-section and branching ratios 
for slepton production, and the number of data and background events kinematically 
compatible for a given mass combination. Signal events have been generated assuming 
model input values of tan/? = 1.5 and /i = -200 GeV/c^. The limits were calculated using 
a likelihood ratio method described in [^. Expected exclusion zones were calculated 
using the same algorithm, from simulated background-only experiments. 

Figure |^.a shows the 95% CL exclusion regions for e^eji production, obtained using 
the full 212 pb~^ of data. For the selectrons, we exclude masses up to mg < 87 GeV/c^, 
providing the mass difference between the selectron and the LSP is above 20 GeV/c^. 
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Figure ||.b shows the 95% CL exclusion regions for P^r^lr production, obtained by 
combining the 183 GeV data with the 189 GeV data. For the smuons, we exclude masses 
up to mr. < 80 GeV/c^, providing the mass difference between the smuon and the LSP 

is above 5 GeV/c^. 

Exclusion limits on ff production were obtained taking into account the signal ef- 
ficiencies for each f — Xi mass point. When determining whether data or background 
events were kinematically compatible with the mass point, the end point of the expected 
momentum spectrum of the visible reconstructed tau was used. Figure ||.c shows the 95% 
CL tr exclusion region obtained by combining the previous data at lower energies with 
the 183 GeV and 189 GeV data, and Figure |[d shows the exclusion regions in the case 
of the mixing angle yielding the minimal cross-section. For the staus, a mass limit can 
be set at 73 to 75 GeV/c^ (depending on mixing) for mass differences between the stau 
and the LSP above 10 GeV/c^. The dedicated search for a low-mass stau yields that a f 
of a mass below 12.5 GeV/c^ is excluded at 95% CL for any mixing angle, provided that 
Am is greater than m^-. 

For the cascade decay analysis, assuming eRCR or p.Rp.R production, one can set ex- 
clusion regions in the SUSY parameter space. To set limits, one has to consider the cross 
section for the selectron or smuon production and the branching ratios for the i 
and X2 ~^ Xil must be taken into account. These cross-sections and branching ratios 
depend on the actual values of the SUSY parameters. The excluded regions for a given 
value of the common scalar mass, mo, and tan/? are presented in Figure( |^) as a function 
of the Higgs superfield mass parameter fi and rhi/2- 

9 Conclusions 

In a data sample of 212 pb~^ collected by the DELPHI detector at centre-of-mass 
energies of 183 GeV and 189 GeV, searches were performed for events with acoplanar 
lepton pairs. The mass limits produced assume input parameters for the Higgs mass 
mixing parameter, fi, of -200 GeV/c^ and ratio of the vacuum expectation value of the 
Higgs doublets, tan/3, of 1.5. 

For the selectron pairs, 56 candidates remained after selection, with an expectation of 
51.2 ±1.5 from SM processes. This allowed a lower limit on the mass for the eR to be 
set at 87 GeV/c^ for Am > 20 GeV/c^. 

In the search for smuon production at 189 GeV, 17 events were selected for regions of 
low Am, with 17.5 ± 0.3 expected from Standard Model processes. For regions of Am > 
35 GeV/c^, 7 candidates were selected, with a background expectation of 9.2 ± 0.2 events. 
At 183 GeV, 5 candidates were selected with a background expectation of 6.1 ± 0.6 
events. Combining these data, a mass limit for fl^ of 80 GeV/c^ was obtained for 
Am > 5 GeV/c^ . 

In the search for stau production, 7 events were selected at 183 GeV with 7.5 ± 0.5 
expected from SM processes. At 189 GeV, 9 candidates passed the selection criteria with 
a background of 10.6 ± 0.7 expected. Combining this data with all our previous data at 
lower energies [0], a mass limit for the stau can be set at 75 GeV/c^ if the stau is purely a 
partner to the right handed tau, and at 73 GeV/c^ if the stau mixing angle is such that 
the production cross-section is minimal. 

In the search for a low-mass stau, 122 events were selected in the mass region below 
15 GeV/c^ and 50 events in the region between 15 GeV/c^ and 27 GeV/c^ . The back- 
ground was 150.1 ± 2.3 and 55.7 ± 1.7 in the two regions. Combining these results with 
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all previous data 0, a f with mass below 12.5 GeV/c^ can be excluded for any mixing 
angle, provided that Am is greater than m,-. 

Events with the topology of e"'"e~77, /i"'"/x~77, and missing energy were analysed and 
a search for C. ix2 ~^ ^IXi "was performed. For the cascade decay selectron search 5 
events were selected with an expectation of 4.2 ± 0.9 from SM process. For the smuon 
case 3 events remained with an expectation of 2.9 ± 0.7. No excess over the Standard 
Model prediction was found. 
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(GeV) 

V V / 


183 


189 


Observed events 


11 


45 


Total background 


12.7 ± 0.8 


38.5 ± 1.3 


Z77 {i^id,ee,TT){n-f) 


1.7 ± 0.2 


1.4 ± 0.1 


4-fermion events 


10.5 ± 0.8 


34.2 ± 1.3 


77 — >-ee, /x/x, rr 


0.5 ± 0.1 


2.9 ± 0.1 



Table 1: Selectron candidates, together with the total number of background events 
expected and the contributions from major background sources. Results shown are for 
54 pb~^ of data analysed at 183 GeVand 158 pb~^ of data analysed at 189 GeV. 





189 GeV 


183 GeV 


Am < 35 GeV/c2 


Am > 35 GeV/c2 


All regions 


Observed events 


17 


7 


5 


Total background 


17.5 ± 0.3 


9.2 ± 0.2 


6.1 ± 0.6 


Z''/7 — > {fj,fj,,ec,TT){nj) 


0.9 ± 0.1 


1.6 ± 0.2 


0.2 ± 0.1 


4-fermion events 


15.7 ± 0.3 


5.5 ± 0.1 


5.8 ± 0.6 


77 — >-ee, /i/j,, tt 


0.9 ± 0.1 


2.1 ± 0.2 


0.1 ± 0.1 



Table 2: Smuon candidates, together with the total number of background events ex- 
pected and the contributions from major background sources. The results are shown for 
the two regions of the SUSY mass space analysed at a centre-of-mass energy of 189 GeV, 
and the full mass spectrum analysed at a centre-of-mass energy of 183 GeV. 



(GeV) 


183 


189 


Observed events 


7 


9 


Total background 


7.5 ± 0.5 


10.6 ± 0.7 


Z77 (/i/i,ee,rr,qq)(n7) 


1.0 ± 0.3 


1.5 ± 0.2 


4-fermion events 


5.0 ± 0.3 


7.2 ± 0.4 


77 — > T'^T" 


0.9 ± 0.2 


0.8 ± 0.3 


77 ^ee, qq 


0.6 ± 0.1 


1.1 ± 0.5 



Table 3: Stau candidates in the search for high mass staus, together with the total 

number of background events expected and the contributions from major background 
sources for centre-of-mass energies of 183 GeVand 189 GeV. 
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TUf < 15 GeV/c^ 


15 GcV/c^ < < 27 GeV/c^ 


(GcV) 


183 


189 


183 


189 


Observed events 


31 


91 


12 


38 


Total background 


38.8 ± 1.5 


111.3 ± 1.8 


14.4 ± 0.9 


41.3 ± 1.4 


ZVt (/^/^,ee,rT,qq)(n7) 


30.9 ± 1.4 


86.5 ± 1.5 


7.6 ± 0.7 


19.9 ± 0.7 


Bhabha 


1.1 ± 0.2 


3.3 ± 0.7 


0.3 ± 0.1 


0.9 ± 0.4 


4-fermion events 


3.9 ± 0.3 


12.8 ± 0.5 


3.9 ± 0.3 


12.6 ± 0.4 


77 T+r-, ee qq 


2.9 ± 0.4 


8.7 ± 0.6 


2.6 ± 0.4 


7.9 ± 1.1 



Table 4: Stau candidates in the search for low-mass staus, together with the total number 
of background events expected and the contributions from major background sources for 
centre-of-mass energies of 183 GeV and 189 GeV. 



Observed events 


5 


Total background 


4.2 ±0.9 


Bhabha 


2.0 ± 0.4 


77 T~^T~ 


1.3 ± 0.7 




0.6 ± 0.3 


4-fermion 


0.3 ± 0.3 



Table 5: Break down of the individual contributions to the tight 6677 sample. Results 
presented are for data taken at centre-of-mass energies of 183 GeV and 189 GeV. 



Observed events 


3 


Total background 


2.9 ±0.7 




0.9 ± 0.4 


e'^e~ — > TT 


1.2 ± 0.5 


4-fermion events 


0.8 ± 0.3 



Table 6: Break down of the individual contributions to the tight /X/X77 sample. Results 
presented are for data taken at centre-of-mass energies of 183 GeV and 189 GeV. 
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e+ ''L,R e+ ^L,R 




^ (a) ll,R e (b) e^,, 

Figure 1: Production diagrams for sleptons in the MSSM. (a) Shows the pair-production 
of sleptons, a possible scenario at collider experiments, (b) Shows the additional t-channel 
contribution to selectron production. 



i £ 




(a) X? (b) X? 

Figure 2: Slepton decay diagrams, (a) Shows the slepton decaying into a lepton of same 
flavour and the LSP. (b) Shows the cascade decay; the slepton decaying into the lepton 
plus the second lightest neutralino, followed by a radiative decay to the LSP. 
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Figure 3: A pre-selection comparison of data and simulated SM events in the selectron 
analysis at 189 GeV. The plots show; (a) Electron pair acoplanarity, (b) Transverse 
momentum of the electron pair, (c) Energy of leading electron, (d) Opening angle between 
lepton pair. The dots with error bars show the data, while the simulation is plotted as a 
histogram. 
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Figure 4: A pre-selection comparison of data and simulated SM events in the smuon anal- 
ysis at 189 GeV. The plots show; (a) Muon pair acoplanarity (b) Transverse momentum 
of the muon pair, (c) Energy of leading muon, (d) Opening angle of the muon pair. The 
dots with error bars show the data, while the simulation is plotted as a histogram. 
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Figure 5: Function used to discriminate against W-pair backgrounds. Left: The points 
represent the data and the sohd histogram shows the contribution from Standard Model 
backgrounds. Right: The discriminating function for a sample of 4-fermion events (solid) 
and SUSY signal events (hashed). Note: the right hand plot is shown for illustration 
purposes only. 
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Figure 6: A pre-selection comparison of data and simulated SM events in the stau analysis 
at 189 GeV. The plots show: (a) Energy of the most energetic, isolated photon (b) The 
square of the transverse momentum with respect to the thrust axis (c) The acoplanarity 
(d) Missing transverse momentum, for events with acoplanarity above 10 degrees. The 
dots with error bars show the data, while the simulation is shown by the histogram. 
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Figure 7: A pre-selection comparison of real data and simulated SM events in the electron 
channel of the cascade decay search. The plots show (a) Momentum of leading electron, 
(b) Energy of leading photon, (c) The invariant mass Mge-^i, (d) The angle between the 
missing momentum and the closest electron. 
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Figure 8: 95% CL exclusion regions for ii production in the MSSM. Figures (a) and 
(b) show respectively the exclusion region for e/j, production in the {i,Xi) mass plane. 
Figures (c) and (d) show the mass exclusion regions for the fji and fmin in the {f,Xi) mass 
plane. The shaded region in the plots shows the obtained exclusion limit, and the solid 
line shows the expected limit treating simulated background as data. In (d), the dotted 
line represents Am = m,-. The limits have been produced using values of tan/?=1.5 and 
= -200 GeV/c^. 
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Figure 9: Exclusion regions at 95 % CL in the SUSY parameter space from the ££77 
events. The shaded region shows the obtained hmit, and the sohd hne shows the hmit 
treating simulated background as data. The exclusion region is obtained assuming a 
slepton mass = 80 GeV/c^, and a value of tan/5 = 1.0. The slepton mass of 80 GeV/c^ 
was chosen as it was the highest excluded mass from the e, fi direct decay search. 




